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Abstract 

Using the OPE, we obtain new sum rules in the heavy quark limit of QCD, 
in addition to those previously formulated. Key elements in their derivation 
are the consideration of the non-forward amplitude, plus the systematic use of 
boundary conditions that ensure that only a finite number of j p intermediate 
states (with their tower of radial excitations) contribute. A study of these 
sum rules shows that it is possible to bound the curvature a 2 = £"(1) of the 
elastic Isgur-Wise function £(w) in terms of its slope p 2 = —£'(1). Besides 
the bound a 2 > |/9 2 , previously demonstrated, we find the better bound a 2 > 
^[4p 2 + 3(p 2 ) 2 ]. We show that the quadratic term |(/3 2 ) 2 has a transparent 
physical interpretation, as it is leading in a non-relativistic expansion in the 
mass of the light quark. At the lowest possible value for the slope p 2 = |, 
both bounds imply the same bound for the curvature, a 2 > j|. We point out 
that these results are consistent with the dispersive bounds, and, furthermore, 
that they strongly reduce the allowed region by the latter for £(io). 
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1 Introduction 



In a recent paper [I] we have set a systematic method to obtain sum rules 
(SR) in the heavy quark limit of QCD, that relate the derivatives of the elastic 
Isgur-Wise (IW) function £(w) to sums over IW functions of excited states. The 
method is based on the Operator Product Expansion (OPE) [2j applied to heavy 
quark transitions [3] and its key element is the consideration, following Uraltsev 
[3], of the non-forward amplitude, i.e. B{vi) — > D^{v') — > B{vf) with in general 
Vi^Vf. Then, the OPE side of the SR contains the elastic IW function £(u>i/) an d 
therefore the SR depend in general on three variables W{, wj and Wif that lie within 
a certain domain. By differentiation relatively to these variables within the domain 
and taking the limit to its boundary, one finds a very general class of SR that have 
interesting consequences on the shape of £(w). 

Let us be more quantitative. As shown in ref. using the OPE - as formulated 
for example in j5] and generalized to Vi ^ Vf 0] -, the trace formalism 6\ and 
arbitrary heavy quark currents 

h = h$T 1 hV , J 2 = fij»r 2 z# (l) 

the following sum rule can be written in the heavy quark limit pQ : 

{ £ ^Tr[B / (^ / )r 2 pW(^]Tr[pWK)ri^(^)]e (n) K)e (n) K) 

D=P,V n 

+ Other excited states} = -2£(w if )Tr [B f {v f )T 2 P[r 1 B i (v l )] . (2) 
In this formula v' is the intermediate meson four-velocity, the projector 

K = \o> + 1>') ( 3 ) 

comes from the residue of the positive energy part of the c-quark propagator, and 
£(wif) is the elastic Isgur-Wise function that appears because one assumes V{ ^ Vf. 
Bi and Bf are the 4x4 matrices of the ground state B or B* meson and T>^ those 
of all possible ground state or excited state D mesons coupled to Bi and Bf through 
the currents. In formula (j2J) we have made explicit the j — ~ D and D* mesons 
and their radial excitations. 
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The variables w i} Wf and w t f are defined as 

Wi = Vi • v' Wf = Vf ■ v' Wif = Vi ■ Vf 

The domain of (w iy Wf, Wif) is [T] 



(4) 



Wi,Wf>l, WiWf - yj(w? - l){wj - 1) < W if < WiWf + yj(w? - l){wj - 1) (5) 
There is a subdomain for Wi = w f = w : 

w>l 1 < w if < 2w 2 - 1 . (6) 

Calling now L(wi,Wf,Wif) the l.h.s. and R(wi,Wf,Wif) the r.h.s. of (J2J), this SR 
writes 

L(wi,w f ,Wif) = R(wi,Wf,w if ) (7) 

where L(wi,Wf,Wif) is the sum over the intermediate -D states and R(wi,Wf,Wif) 
is the OPE side. Within the domain (0) one can derive relatively to any of the 
variables Wi, Wf and Wif 

Qp+q+rl d p+q+r R 

dwfdwjdwif dwidwjdwif 

and obtain different SR taking different limits to the frontiers of the domain. One 
must take care in taking these limits, as we point out below. 

Let us parametrize the elastic Isgur-Wise function £(w) near zero recoil, 

aw) = l-p 2 (w-l) + ^(w-l) 2 ---- (9) 

From the SR (j2j), we gave in ref. pQ a simple and straightforward demonstration of 
both Bjorken [7] jS] and Uraltsev jlj SR. Both SR imply the lower bound on the 
elastic slope 

p 2 = -£'(i) > ~ • (io) 

A crucial simplifying feature of the calculation was to consider, for the currents 
(Jljl. vector or axial currents aligned along the initial and final velocities i;, and Vf. In 
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we did also obtain, modulo a very natural phenomenological hypothesis, a new 
bound on the curvature. 

- 2 = ni)>^ 2 >^- (ii) 

This bound was obtained from the consideration in the SR of the whole tower of j p 
intermediate states [§]. A crucial feature of the calculation was the needed derivation 
of the projector on the polarization tensors of particles of arbitrary integer spin [TU] . 

Using the SR involving the whole sum over all j p intermediate states, we pursued 
our study in jTT] and did demonstrate that the IW function £(w) is an alternate series 
in powers of (w — 1). Moreover, we did obtain the bound for the n-th derivative at 
zero recoil (-l) n £^(l) 

i-ir&Ki) > ^±Vir^) > (12) 

demonstrating rigorously the bound (fTTjl and generalizing (JTU|) and (|TT|) to any 
derivative. 

The aim of this paper is to investigate whether a systematic use of the sum rules 
can allow to obtain better bounds on the curvature. As we will see below, the answer 
is positive. The reason is that only a finite number of j p states, with their radial 
excitations, contribute to the relevant sum rules and one is left with a relatively 
simple set of algebraic linear equations. As we will see, this is due to the crucial 
fact that we adopt particular conditions at the boundary of the domain (JSJ). 

2 Vector and Axial Sum Rules 

We choose as initial and final states the B meson, 

Biivi) = P J+ (- 75 ) Bf{v f ) = P/+(-T5) (13) 

where the projectors P; + , Pf + are defined like in Q. Moreover, we consider vector 
or axial currents are projected along the Vi and Vf four- velocities. Choosing the 
vector currents 
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Ji = h ( :h t h^ , J 2 = T$j1> f h!$ (14) 

and gathering the formulas (48) and (89)-(91) of ref. PP we obtain for the SR (J2J) 
with the sum of all excited states j p , as written down in JT] : 



w, 



l>0 n 



+ E S e( W i> W f' W if) E r i-S( W r i-l/2( lW /) = (1 + + U7/ + w if )£(w if ) . (15) 



1/2 V" 7 '*; '£-1/2 
€>1 n 



Choosing instead the axial currents 



Ji = t&^h$ , J^hSjtj-fth® (16) 
the SR (j2J) writes, from the formulas (48) and (92)-(94) of ref. pQ, obtained in JT] : 



£>0 

-(1 - - W/ + Wif)Z[w if ) . (17) 



+ (to< - l){w f - 1) E w /s Wi/) E r i-i ( /2K) r i-i ( /2( w /) 



Following the formulation of heavy-light states for arbitrary j p given by Falk [9 , 
we have defined in ref. [T] the IW functions T^y\ (w) and rf}^\{w)^ that correspond 
to the orbital angular momentum i of the light quark relative to the heavy quark, 
j — I ± ■= being the total angular momentum of the light cloud. For the lower values 
of £, one has the identities with the traditional notation of Isgur and Wise jH] : 



r[%(w) = aw) , t[ 1 / 1(w) = 2t 1/2 {w) , rff 2 (w) = V3 r 3/2 (w) (18) 

(i) 

1/2 



where a radial quantum number is implicit. Therefore, the functions t[]\(w) and 



Ty\(w) correspond, respectively, to the functions C(w) and t(w) defined by Lei- 
bovich, Ligeti, Steward and Wise [12]. 

In equations (fT3|) and (fTT)|) the quantity S n is defined by 

= v fui ■ --v fUn T"' «^ • • -v ifln (19) 
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and the polarization projector r p yi "' Uk, ^ 1 "'^ n , given by 



rpu 1 ---U n ,fl 1 ---/J, n _ £ /(A)*fl —v„ £ t(X)ni---ti n (20) 

A 

depends only on the four-velocity v' . The tensor e'Wm-Vn j s ^he polarization tensor 
of a particle of integer spin J = n, symmetric, traceless, i.e. £ / ^ M1 "" Atn <7^- = 
(i ^ j < n), and transverse to v', v' e'^^ 11 '"^ = (i < n) [I] [TQJ . 

Moreover, as demonstrated in the Appendix A of ref. [I], S n is given by the 
following expression : 



S n (wi,w f ,Wi f ) = J2 C nik {wl-l)\w)-l) k {w i w f -w if ) n - 2k (21) 

with 



0<fc<§ 



k (n\) 2 (2n-2k)\ 



Cn > k ~( 1} (2n)! fc!(7i-Jfe)!(n-2Jfe)! ' (22) 



The relation 



L (wi,w f ,Wif) = R (wi,w f ,Wif) (23) 

Wif = l,Wi=Wf=W Wif = l,Wi=Wf=W 

gives, dividing by 2(w + 1), Bjorken SR [7] jH], including now the whole sum of 
intermediate states : 



£>0 
W — 1 



2 ^>o ^ + 



+ VE ^ - I)*" 1 E ^Hrj^H = 1 (24) 



where, from (J2U) and (|22jl 



W) 2 



S n (w,w,l) = C n (w 2 -l) n C n = C n , k = 2 n ^-. (25) 

o<fc<f ^ Zn >- 

Remember that, usually, the first terms in the sum (J24|) are written in the notation 
(ITSIl of Isgur and Wise : 



tu + 1 



(«), 



+ (w - 1 



E 



r (n) 
r l/2 



(w) 



+ (w + l) 



T. 



in). 
3/2 1 



W 



+ 



1 . 

(26) 
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Going now to the axial current SR (fT7|). the condition 



L A (wi,w f ,w if ) = R A (wi,w f ,w if ) (27) 

" J 111. .. 1 Til- 111 ,. Til V " J ' 1(1. n 1 111. 111 111 V ' 



Wif =1 ,Wi=Wf =w 



Wif=l,Wi=W f=W 



gives again, dividing this time by 2(w — 1), the complete Bjorken SR (j23)- Notice 
that, as it must, one obtains the same SR from the vector (J 14)) and the axial current 
(ITTIj) because, from (|23j) . one has 

(2n+l)C n+1 = (n + l)C„. (28) 

3 Equations from the Vector Sum Rule 

In what follows, to look for independent relations, we make use of the fact that 
the SR (fTo'j) and (|T7j) are symmetric in the exchange Wi <-» Wf. 

Let us first consider the derivatives of the SR for vector currents (115)) relatively 
to Wif with the boundary condition Wif = 1. For Wif = 1, the domain © implies : 

Wi = Wf = w . (29) 

We define therefore 

L v (w if ,w) = L v (wif,w h Wf)\ Wi=wf=w 

R v (w if ,w)=R v (w if ,Wi,w f )\ Wi=Wf=w . (30) 

We then take the p + q derivatives 

( dP + "L v \ ( dP+Wv \ rqi , 

(31J 



\ dwf f dw q / \ dwfrdwi j 

and exploit systematically the obtained relations. To get information on the curva- 
ture a 2 of the elastic IW function Q we need to go up to the second order derivatives. 
As we will see, a crucial feature of the adopted boundary condition is that only a 
finite number of j p states, with their tower of radial excitations, contribute to the 
sum rule. Notice that we could have derived first relatively to Wi and take the limit 



7 



Wi = 1, and then derive with respect to w = Wif = Wf. We do not obtain however 
new information from these sum rules that with the former boundary conditions. 

Let us proceed with care and begin with the first order derivatives. From (JTSJl 
and (}3T]) . we obtain the following results. 

With the notation (jSTJ), for p = q = we obtain the trivial result £(1) = £(1), 
while for the derivatives p = 1, q = we obtain Bjorken SR for the slope p 2 : 



(32) 



A ' q ^ ^'J 1 4 ^ V V 2 

° n n 

The relation to the Isgur-Wise notation is given by (|18|) . 

Going now to the second order derivatives, we find the following relations. For 
p = 0, q = 1 we get £(1) = £(1). For a purpose that will appear clear below, we 



make explicit the IW functions between j p 

and Wise £ (n) H <P1>- 
For p = 2, q = : 



states using the notation of Isgur 



1o y 
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p= 1, q= 1 : 



(2)(n) 
5/2 



(1) 



+ 



(33) 



)(n). 



3/2 



~E^ 

U n 



(!)(*)/ 



3/2 (^) T 3/2 



(!)(«)'/ 



(i)-E^ (B) (i)^ w '(i 



(!)(«)', 



- 2 E[ 



r$ w (i; 



24 



.(2)(n) 
5/2 



(1) 



(34) 



p = 0, g = 2 : 



1 - 8p 2 + 4a 2 + 4 £ [£W (1)] 2 + 8 E [t$ (b) (1)] " + E [t$ (b) ( X 



+ ?E 



o ^ T 3/2 \ L ) T 3/2 { 



i)+4S><;> l ^i)Tm ( ">'(i) 



(2)(n) ( 



3 ^L r 3/2 '(1. 



.(2)(n) 
5/2 



(1) 



(35) 



8 



The equations (J33J)-(J33]) are a set of linear equations in the elastic slope p 2 and the 
curvature a 2 , and the following quantities, that are series on the radial excitations, 
indicated by the sums over n : 



n 



.12 



3/2 1 



-E 



3/2 



(1)t, 



3/2 1 



n 



.(2)(n). 
3/2 1 

.(2)(n) 
5/2 1 



(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 



We realize first that, due to the fact that we compute the second derivatives in (}3*Tj) 
{p + q = 2) and use the boundary conditions Wif = w = 1, the series in j p states is 
truncated and includes at most the £ = 2 states j p = | , | , corresponding to the 
unknowns (|41|) and (|4*2*|) . On the other hand (|36|) is the square of the derivatives at 
zero recoil of the lowest j p = ~ , and (|37|) and (|38|) depend on the IW functions 
of the transitions to the P-wave states j p = | + , | + , that are simply related to the 
slope p 2 through Bjorken and Uraltsev SR, as we write down below again. Finally, 
we have two other unknowns (J3*9*|) and (|40|) that involve the derivatives of the P- 
wave IW functions r^l (w), Ty 2 (w) at zero recoil. These quantities were already 
introduced in ref. pQ. 



4 Equations from the Axial Sum Rule 

Let us now consider likewise the derivatives of the SR for axial currents (j!7|) 
with the boundary condition w^f — 1, Wi — Wf—w— > 1 : 



d p+q L A \ ( d p+q R A 



dw P f dwi j \ dw P f dw q j 

>■ l J ' W i f=W = l V l J ' Wif =w = l 



(43) 
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Since all terms in (j!7|) vanish for W{ = Wf = Wif = 1, to obtain information on the 
curvature a 2 we will need to go up to the third order derivatives. 

With the notation for p = q = we get = 0, and for p — 1, q — and 
p = Q,q=l, £(1) = £(1). 

For the second order derivatives we obtain the following results. For p = 2, g = 
and p = g = 1 we get the same relation 



while for p = 0, q = 2, we get Bjorken SR : 



(44) 



Both equations ({H| and (f43|) imply Uraltsev SR 



2 1 



(45) 



2 1 

~ 4 



[ r l/2 



2 _ 1 

~ 4 



(46) 



using the notation (fTH|) . 

Going now to the third order derivatives, we obtain the following results. 
For p = 3, q = : 



2EhT"(i) 



p = 2, g = 1 : 



^ = 2E^ B) (l)^ B) '(l)+6E[r^ (n) (l) 



p= 1, ? = 2 : 



(47) 



(4£ 



^E[^w] 2 +2EB w (i) 



+ BErg W (l)^ W (l) + ^EK W (1)] 2 + f £ [r$ W (l) 



(49) 



p = 0, q = 3 : 
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-3p 2 + 3a 2 + 3E[e ( " ) '(l)] 2 + 4E[r 3 ( / i r ) (l) 



+ 8E 



.(!)(")■ 
3/2 1 



(!)(«)' 



(1) 



+ 2E[^ w (i)] a + ^E[^ (B) (i 



. 



(50) 



The equations (|37)l - (j5U|l depend on p 2 , a 2 and the same set of unknowns listed in 
Pfi|) - (j4"2|l . Let us now look for the set of linearly independent equations that can be 
obtained from the sets of equations obtained from the Vector and Axial SR. 



5 Linearly independent relations 

Let us concentrate on the equations (}33|) - (j35|) and (j4"T|) - (j50j) obtained respectively 
from the Vector and Axial Sum Rules. 
Using Bjorken SR 



(61) 



the relation 



-J>~ 



i = EK (n) W 



(52) 



obtained from (|32j) and (|44j) . and f|44|) and (J47)) we obtain finally the following set 
of relations that are linearly independent 



p 2 = - | W(i)^ w (i) 



(53) 



a 



^ = 2£[r 5 ™»>(l) 



(54) 
(55) 
(56) 
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^ 2 -^ 2 +E[^'(i) 



i2 



. 



(57) 



3 3 

Relations (jo^j) and (p>fj) were obtained in ref. [I], and relation (|B3j) and (JoT)j) in 
[TTj . The systematic study of the present paper using all possibilities (pTTj) and 
involving the curvature gives the new equation (j!T7j) . 



6 Bounds on the curvature 

The last two equations ()56|) and (|57p involve the curvature with a negative sign 
and positive definite quantities. Making explicit in the sum £[£ {n) '(l)] 2 the ground 

n 

state IW function slope (1) = — p 2 , one obtains the two equations 

p 2 -^ 2 + EK ( /1 (n) (l)! 2 = 0. (58) 

n 

|p 2 + (p 2 ) 2 -^ 2 + Eie (n) '(l)| 2 = (59) 
that imply respectively the bounds : 

* 2 > \ P 2 (60) 



2 1 

a 2 > - 



V + KP 1 ) 2 • (61) 
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The bound (|5U|) was obtained in ref. [T| using the relations and (J53j) and making 
the assumption 



-E^ (B) (1)^ {B)/ (1)>0- (62) 

n 

Later, (J6Uj) was demonstrated rigorously in ref. ^T] and generalized to the n-th 
derivative. However, in this latter paper, only derivatives relatively to Wif were 
taken, while in the present work a systematic use of (J3T|) and (pfHj) is carried out. 

The inequality (JSTjl is the best of the bounds that we have obtained for a 2 for 
any value of p 2 , and is the main result of this paper. 

Interestingly enough, both bounds (fHUJ) and (jfilj) coincide at the lower bound 
p 2 > | implied by Bjorken and Uraltsev SR, and (jjHj) . At the value p 2 = \ one 
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then gets indeed the same absolute bound (i.e., independent of p 2 ) for a 2 , namely 



7 Implication on the P-wave IW functions at zero 
recoil 

Let us now express the sums of products of the P-wave Isgur-Wise functions 
i- _> 1+ an d their derivatives ErS (n) (l)rg (n) '(l) and £T$ (n) (l)r$ (n) '(l) in 

n n 

terms of p 2 and a 2 . From (|5Hjl and (jf>4j) we obtain, using now the notation of Isgur 
and Wise (HHJ) : 

-E^(l)<(l) = k (63) 

-E-i ( ; 2 ) aK ( /t(i) = -^p 2 +^ 2 - (64) 

Using the bounds (jlUj) and (fTT|) for p 2 and cr 2 one finds 



-E^(l)^'(l)>4 (65) 
-£t$(1)t$(1)>0. (66) 

n 

Strictly speaking, these relations do not give information on the slope of the 

fnV fnV 

lowest n = IW functions (1) and (!)• However, if the n = state dominates 
the sum, the inequalities (|63j) and (|6^j) imply that the slopes t^9 2 (1) and rfl (1) are 
negative, as it is plausible on physical grounds for form factors that do not involve 
radially excited states. 

This is indeed the case for the Bakamjian-Thomas type of quark models, that 
satisfy IW scaling [T3] and Bjorken and Uraltsev sum rules [Hj. We have conjectured 
in pQ that this class of models presumably satisfy all the SR of the heavy quark limit 
of QCD. 

In the Bakamjian-Thomas model one finds for the phenomenologically successful 
spectroscopic model of Godfrey and Isgur [T^], the numbers 
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- (1) = 0.43 (67) 

-r^(l)rf^(l) = 0.04 (68) 

that by themselves satisfy the preceding bounds, so that the n = state seems to 
give a dominant contribution to the l.h.s. of (ffiHJ) and (jHEJ)- 

8 Non-relativistic limit of the bounds 

There is a simple intuitive argument to understand the term |(p 2 ) 2 in the best 
bound (}6T| . Let us consider the non-relativistic quark model, i.e. a non-relativistic 
light quark q interacting with a heavy quark Q through a potential. The form factor 
- to be identified with the IW function - has then the simple form : 

F(k 2 ) = / dr ^(r)exp(i — ^ — k-r)^ (r) (69) 
J \ m q + rriQ J 

where <fo{ r ) is the ground state radial wave function. In the small momentum 

transfer limit, the IW variable w writes, in the initial heavy hadron rest frame : 

v' 2 k 2 

Identifying the non-relativistic IW function £,nr{w) with the form factor F(k 2 ) (JBTJj) . 
one finds, because of rotational invar iance : 



£ NR (w) = l-m 2 q < 0\z 2 \0 > (w - 1) + - - m\ < 0\z 4 \0 > (w - l) 2 + • • • (71) 

where |0 > stands for the ground state wave function, and we have neglected in 
the (w — l) 2 coefficient subleading terms in powers of l/(m q z) (internal velocity). 
Therefore, one has the following expressions for the slope and the curvature, in the 
non-relativistic limit : 



p 2 NR = m 2 q <0\z 2 \0> a 2 NR = ~ m\ < 0\z 4 \0 > . (72) 



From spherical symmetry one has 
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< 0|^ 4 |0 > = - < 0|r 4 |0 > 
5 



Using now completeness \ n >< n \ — 1; 



< 0|r 4 |0 > = | < 0|r 2 |0 > | 2 + ^ I < ^|r- 2 |0 > | 2 



we use again spherical symmetry 



(73) 



(74) 



<0|r 4 |0> =9| <0|2 2 |0> | 2 + 9 \<n\z 2 \0>\ 2 

nj^0,rad 

where the latter sum runs only over radial excitations. 

Therefore, from (j72|) - (f73j) we can rewrite (J 2 NR under the form 



(75) 



or 



a 



NR 



and therefore 



3 
5 



m 2 | < 0|z 2 |0 > 



+ m q I < ™k 2 |o > 

n^0,rad 



a 



NR 



2 

Pnr 



2 3 4 v- 
+ 5 m i L 

nj^0,rad 



< n\z 2 \0 > I 2 



(76) 



(77) 



a NR > - 



3 

5 



2 

Pnr 



(78) 



Notice that, denoting by R the bound state radius and m q the light quark mass, 
in the non-relativistic limit, just from expressions (J72j) . one can see that p 2 NR scales 
like m 2 R 2 , while a 2 NR scales like m 4 i? 4 and both the l.h.s. and the r.h.s. in (|78j) 
scale in the same way. 

Going back to the relativistic bounds ()6U| ) -([6i p . we observe that the terms pro- 
portional to p 2 are subleading in the non-relativistic expansion and correspond to rel- 
ativistic corrections specific to QCD in the heavy quark limit. In the non-relativistic 



limit p 2 



~ m 



2 R 2 3> 1, and the power (p 2 ) 2 is leading. We can understand therefore 



the appearance of the term l(p 2 ) 2 in the r.h.s. of the inequality (pTTj) . 
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9 An example of fit to the data 



An interesting phenomenological remark is that the simple parametrization for 
the IW function [T3| 

/ 2 \ 2p2 

«<•> - (^) p»> 

gives 

* 2 = P Y + (p 2 ) 2 (so) 

that satisfies the inequalities (jS^-flU) if p 2 > §, i.e. for all values allowed for 
p 2 . Moreover, interestingly, at the lowest bound of the slope p 2 = |, (fKU|) implies 
precisely the lower bound of the curvature a 2 = y|, as pointed out in [TT] . 

Notice that, in ref. [T^], within the class of Bakamjian-Thomas quark models, 
the approximate form (J79j) was found with p 2 = 1.02 in the particular case of the 
spectroscopic model of Godfrey and Isgur. This gives a curvature ()8()|) a 2 = 1.55, 
close to the bound (JoTj) . that gives a 2 > 1.44, stronger than the bound (JtjUj) . that 
implies cr 2 > 1.27. 

As a simple example of a fit with the simple function (fTH|) . we can use BELLE 
data on B° — > D* + e~v for the product | V c b\J r *{w) [T2], as shown in the Figure. 

The function J r *{w) is equal to the Isgur- Wise function £(w) in the heavy quark 
limit. Assuming only departures of this limit at w = 1, i.e. fitting £,(w) from the 
data with 

\V ch \T*{w) = \V db \F'(l)S(w) (81) 
we obtain the following results for the normalization and the slope : 

.F*(l)|y c6 | = 0.036 ±0.002 p 2 = 1.15 ±0.18 (82) 

with the other derivatives of £(w) fixed by (|T§j) (Fig. 1). 

As we can see, the determination of ^-"*(1)|\46| is rather precise, while p 2 has a 
larger error. However, the values obtained for | V c b\ and p 2 are strongly correlated. It 
is important to point out that the most precise data points are the ones at large w, 
so that higher derivatives contribute importantly in this region. Due to the alternate 
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Figure 1 : Fit to T* (w) \ V& \ using the phenomenological formula f7D|) and the BELLE data for 
B — > D*iv [16], assuming only violations to the heavy quark limit at w = 1. The fit gives the 
results i)8'2l) . 
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character of £(«?) as a series of (w — 1), one does not see strongly the curvature of 
£(w) in the Figure, but the curve is definitely not close to a straight line. Linear 
fits, as are commonly used, should be ruled out at the view of the bounds that we 
have found. 

We must emphasize that the fit that we present is a simple exercise in the heavy 
quark limit. Radiative corrrections and I/ttlq corrections, that enter in the relation 
between the actual function F(w) and its heavy quark limit should be taken 

into account, although this does not seem to be an easy task |17j . 

10 Comparison with the dispersive bounds 

A considerable effort has been developed to formulate dispersive constraints 
on the shape of the form factors in B — > D*lv [IHj-j221- The starting point are 
the analyticity properties of two-point functions and positivity of the corresponding 
spectral functions. Then, dispersion relations relate the hadronic spectral functions 
to the QCD two-point functions in the deep Euclidean region, and positivity allows 
to bounds the contribution of the relevant states, leading to constraints on the 
semileptonic form factors. 

We will now compare our method, that gives information on the derivatives of 
the Isgur-Wise function, with the dispersive approach. 

A first remark to be made is that our approach, based on Bjorken-like SR, holds 
in the physical region of the semileptonic decays B — > D^*'iv and in the heavy quark 
limit. However, concerning this last simplifying feature, we should underline that 
there is no objection of principle to include in the calculation radiative corrections 
and subleading corrections in powers of 1/rriQ. 

The dispersive approach starts from bounds in the crossed channelby comparison 
of the OPE and the sum over hadrons coupled to the corresponding current, BD, 
BD* , ■ ■ ■ Then, one analytically continues to the physical region of the semileptonic 
decays. This is done for a single reference form factor, for example the combination 



that enters in the B — ► Div rate. In the heavy quark limit h-(w) = 0, V\(w) = 




(83) 
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h + (w) = £(w). Ratios of the remaining form factors to V\(w) are computed in the 
physical region by introducing 1/mq and a s corrections to the heavy quark limit. 
The dispersive approach considers physical quark masses, in contrast with the heavy 
quark limit of our method. 

The two approaches are quite different in spirit and in their results. However, it 
can be interesting to compare numerically our bounds with the ones of the dispersive 
approach, as they happen to be complementary. We must however keep in mind 
precisely the differences between the two methods. 

We have demonstrated in ^T] that the IW function £(u>) is an alternating series 
in powers of (w — 1), with the moduli of the derivatives satisfying the bounds (II 2j) 
and JBTJ). 

10.1 Comparison with the work of Caprini, Lellouch and 
Neubert 

Let us consider the main results of ref . j2I| , that are summarized by the one- 
parameter formula 

Y^l 1 - 8p 2 z + (hip 2 - 10> 2 - (252p 2 - 84> 3 (84) 
Vi(l) 

with the variable z(w) defined by 



z = - = (85) 



and the allowed range for p being 



- 0.17 < p l < 1.51 . (86) 
Of course, the function Yj ^) contains finite mass corrections that are absent at 

Vi(l) 

present in our method. Nevertheless, let us first compare these results with our 
lower bounds (|12|). assuming the rough approximation 

1H W - (87) 

Of course, since the expansion (|Mj) stops at third order in z, it would only make 
sense in the comparison to go up to the third derivative of £(w). Using our notation, 
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the results of Section 4 of ref. [21] for the first derivatives write, from the expansion 
(JEIJ), in terms of the slope p 2 : 



£"(1) = l(67p 2 - 10) 
ni) = -^( 148 V- 372) 



(89) 



8. 



>8 



with p 2 in the range (jEEJ). From (j88|) and (J89|) . using the notation £(w) = 1 — p 2 {w — 
l)+c(w — \) 2 +d{w — 1) 3 + - • • one gets the numerical relations |2I] : c = 1.05p 2 — 0.15, 
d = -0.97p 2 + 0.24. 

The lower bound (|86|) on p 2 is very weak : there is a region for p 2 that is below 
the Bjorken bound p 2 > 1/4, and this reflects on the regions allowed for the higher 
derivatives. Of course, we must always keep in mind that in these considerations we 
are neglecting finite mass and a s corrections. Moreover, one should notice that the 
series is not alternate in the whole range ()86|) . We have verified analytically that 
the expression (jElJ) is an alternate series in powers of (w — 1) for values of p 2 > |, 
i.e. for p 2 satisfying precisely Bjorken bound. It is not clear to us whether this is 
just a numerical coincidence or it has a deeper significance, i.e. if it extends to the 
complete series in powers of z. 

Let us now comment on the implications of our bounds (fT2*J) . The first important 
remark is that, within the simplifying hypothesis f)87j) . the range ()8(j)l is considerably 
tightened by the lower bound on p 2 > | implied by Bjorken and Uraltsev sum rules. 
Therefore, we will consider hereafter, instead of (j86|) . the improved range 



that shows that our type of lower bounds are complementary to the upper bounds 
obtained from dispersive methods. Within the hypothesis of the heavy quark limit, 
the region allowed by the dispersive bounds for £(w) with p 2 within the range (jHHJ) 
is obviously much reduced by the bounds (|9T)j) (Fig. 2). 

On the other hand, it follows that the second and third derivatives (J85j) . (jBHJ) 
satisfy the bounds (fT^j) 



3 9 

- < p 2 < 1.51 
4 



(90) 



= ^(67p 2 - 10 > \p 



(91) 
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Fi gure 2: Tiie upper (lower) curves are the representations of t;(w) according to the dispersive 
approach of Caprini et al [21] (84 )) -i|#7 )) . The upper (lower) curve correspond to p 2 = —0.17 
(p 2 — 1.51). The shadowed region is the region forbidden by the Uraltsev bound p 2 > |. The 
remaining allowed region corresponds to The curve within this allowed region is the best fit 

to BELLE data [16], normalized to w = 1, that gives T* (l)\V cb \ = 0.036 ±0.002, p 2 = 1.16 ±0.15, 
in practice the same fit as \8'2fy with the phenomenological formula 



-r (1) = ^( 14 8V - 372) > ^"(1) = I ^(67p 2 - 10) (92) 

respectively for p 2 > 0.36 and p 2 > 0.42. Therefore, for values of p 2 that are within 
the range (|9"D|). these inequalities are a fortiori satisfied, but they are not satisfied 
in the whole range (|86|). 

Finally, let us look for the implications of our improved bound on the curvature, 
eq. (|6ip. Combining the linear dependence obtained from dispersive methods (|88|) 
with the inequality (jol*|) one obtains the condition 

^(67p 2 -10)>^[4p 2 + 3(p 2 ) 2 ] (93) 

that gives the range 
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0.28 < p 2 < 1.88 . 



(94) 



Interestingly, the condition (|93|) gives by itself an upper bound for p 2 that is of 
the same order than the upper bound (|86|). Moreover, the range ()94|) contains the 
improved range (J90|) . and things appear to be coherent. 

10.2 Comparison with the work of Boyd, Grinstein and 



Let us now compare with the dispersive method results of the work of Boyd, 
Grinstein and Lebed [221 • in this work, the QCD part of the calculation includes a s 
and non-perturbative (condensate) corrections, and new poles below the annihilation 
threshold ignored in [21]. In a form that allows to make the comparison with our 
results, the authors of [22] obtain the following expansion for the scalar form factor 



Lebed 



/o(l) + [1.72a! - 0.77/„(l)J (w - 1) 

+ [-1.74a! + 0.21a 2 + 0.55^,(1)1 (w - l) 2 + • • • 



(95) 



where fo(w) has been defined by Caprini and Neubert j2HI 



fo(w(q 2 )) 



fo(q 2 ) 



(96) 



(M B - M D )y/M B M D (w + 1) 



with 




(97) 



f±(q 2 ) being the form factors governing the rate B Dtv (q = p — p 1 ) : 



< D(p')\V,\B(p) >= f+{q 2 ){p + p')» + f~{q 2 ){p-p')» . 



(98) 



Heavy quark symmetry implies 



Uw{q 2 )) = i{w) . 



(99) 
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The coefficients a n in are defined by the expression for a generic form factor 

ini, m 



-1 oo 

F{z ^rnm^ a " z " (100) 

where z is defined by ()85|) . The functions -P(z) and (^(2;) - respectively the Blaschke 
factor and the outer function - contain the subthreshold singularities in the annihi- 
lation channel, respectively the B c poles and the kinematic singularities. The basic 
result of the dispersive approach is that the coefficients a n of the series obey 

00 

X>n<l- (Id) 

n=0 

To compare with our bounds we proceed like we did above. Since in the heavy 
quark limit (}99|) holds, we set /o(l) = 1 and write the Isgur-Wise function in terms 
of the coefficients a„ : 



£(w) = 1 + (1.72a! -0.77)(w- 1) 

+ (-1.74a! + 0.21a 2 + 0.55) (w - l) 2 + • • • (102) 

Notice that in (jHHj) and (J102)) it does not make much sense to consider higher powers 
(w — l) n (n > 3) unless the corresponding a n (n > 3) are introduced. Then, our 
lower bounds (|T2*|) write 

-1.72O1 + 0.77 > 1 

2 (-1.74ox + 0.21a 2 + 0.55) > -(-1.72ai + 0.77) (103) 
implying, respectively 



ai < 0.01 

a 2 > 3.17ai — 0.33 



(104) 



Since, from ()104j) and (jlOlj) we have 
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- 1 < oi < 0.01 (105) 
and the coefficient of ai in (|104|) is large, the whole range 

- 1 < a 2 < 1 (106) 

is allowed. This seems to support the statement of ref. [22] that a 2 cannot always 
be neglected. 

Moreover, using now the quadratic bound (|6ip. one obtains 

3(p 2 ) 2 -6p 2 + 2(l-a 2 ) < (107) 

and therefore 

-0.5<a 2 <l (108) 
giving the following range for p 2 in terms of a 2 : 



!-^<^< 1 + ^ (109) 

giving the wide range 

< P 2 < 2 . (110) 
For a 2 = 0, implicitely assumed in ref. [21], one finds the range 



0.42 <p z < 1.58 (111) 

a domain qualitatively consistent but somewhat narrower than the corresponding 
one (|94|) obtained from the linear relation between the curvature and the slope given 
by ref. [21]. 

In conclusion, there is no contradiction between the dispersive bounds and the 
type of bounds that we have obtained using Bjorken-like sum rules in the heavy 
quark limit. The latter appear rather as lower bounds that are complementary to 
the upper bounds of the dispersive approach, tightening considerably the allowed 
range for p 2 and for the higher derivatives of £,{w) as well. 
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11 A phenomenological Ansatz for the Isgur-Wise 
function and the dispersive constraints 



At the light of the preceding discussion, we are now going to address the question 
of whether the phenomenological Ansatz for the IW function proposed in Section 9, 




(112) 



satisfies, assuming the heavy quark limit (J%7|) or (ffijj) . the constraints of the disper- 
sive approach. 

We will follow here the formulation of Boyd et al. [22J and consider the form 
factors f+(q 2 ) and fo(q 2 ) defined by (|H7jl and |J9B}) . In the heavy quark limit, one 
has the relations 



u( q \w)) - M f + Md e(«o (us) 

f (q 2 (w)) * (M B - M D )y/M B M D (w + . (114) 

We denote now generically any of these form factors by F(q 2 (w)), or through 
the transformation (jHSJ), F{q 2 {z))- 

We adopt the phenomenological formula ()112|) for £(u>) and define the corre- 
sponding series ()100|) 



oo 



a n z n = P(z) <p(z) F(z) (115) 

n=0 

where P(z) and (p(z) are the Blaschke factor and the outer function of the corre- 
sponding form factors. 

We want now to compare the coefficients a n obtained from (jll3j) - (jll4|) . assuming 
F(z) = £(w(z)) given by (fTT2l . to the condition (ITTTTJ) 

oo 
n=0 

The outer functions (p(z) and the Blaschke factors P(z) for f+(q 2 ) and fo(q 2 ) are 
given in ref. j22j, respectively by formula (4.23) and Table 1 and by formula (4.25) 
and Table 3. We have singled out f+{q 2 ) and fo(q 2 ) as given by ()113|) and ()114|) but 
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we could have taken any other form factor related, up to kinematic factors, to £(w). 
Of course, the results for the coefficients a n would differ according to the considered 
form factor. 

We use the numerical parameters of this paper, and two choices for p 2 in formula 
(I112j) . namely p 2 = 1.023, that corresponds to the Isgur-Wise function obtained 
within the Bakamjian-Thomas scheme from the Godfrey-Isgur spectroscopic model, 
as found in ref. |T3], and p 2 = 1.15 obtained from the fit of Section 9. 

Denoting the Blaschke factor and outer function for each form factor by the 
corresponding subindices, we find, for p 2 = 1.023, the series ()115j) for f + (q 2 (z)) : 

P + (z)(p + (z)f + (q 2 (z)) = 0.0143 - 0.0179,2 -0.1164/ + 0.3277/ 

- 0.1995/ - 0.4497/ + 1.2347/ + • • ■ (117) 

and for f Q (q 2 (z)) : 

P (z)(po(z)fo(q 2 (z)) = 0.0834 -0.1750,2 - 0.1725/ + 0.8673/ 

- 1.1600/ - 0.8943/ - 0.4346/ + • • • (118) 

For p 2 = 1.15 we find, respectively : 

P + (z)(p + (z)f + (q 2 (z)) = 0.0143 - 0.0326,2 - 0.0907/ + 0.4294/ - 

- 0.5779/ - 0.0306/ + 1.3868/ + • • • 

and 

P (z)ipo(z)f (q 2 (z)) = 0.0834 - 0.2599,2 + 0.0484/ + 0.9094/ - 

- 2.0079/ + 2.5089/ - 2.3596/ + - • • 

Comparing to the condition ()116|) . we observe two points. First, the first coeffi- 
cients have squares well below 1, specially for f + (q 2 ). That this happens to be the 
case for this form factor, that has three Blaschke factors, reinforces the idea that 
one should be closer to the IW function, according to our hypothesis (H13J) - (jll4J) . as 
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the number of Blaschke factors increases. Second, high powers of z have coefficients 
that can be of 0(1) and they are strongly dependent on the value of p 2 , specially 
for fo(q 2 ), that has only two Blaschke factors. We have actually observed that the 
behaviour of the coefficients oscillate, as can be seen in ([118)1 . 

Our conclusion is that, owing to the fact that the coefficients, up to order z 3 
included, satisfy the condition fflKij) . the "dipole" formula ()112j) gives, on phe- 
nomenological grounds, a good enough representation of the form factors (|113|) . 

(hud. 



12 Conclusions 

In conclusion, using sum rules in the heavy quark limit of QCD, as formulated 
in ref. pQ, we have found an improved bound for the curvature of the Isgur-Wise 
function o 2 = £"(1) > §[4p 2 + 3(p 2 ) 2 ] that implies the already demonstrated [I] [TT] 
absolute bound a 2 > y|. 

Beyond the simple Ansatz (|79|) introduced above, any phenomenological parame- 
trization of £(u>) intending to fit the CKM matrix element | 1 in i5 — ^ D^*Hu should 
have, for a given slope p 2 satisfying the bound (fTUj) . a curvature a 2 satisfying the 
new bound (JEH). 

Moreover, we discuss these bounds in comparison with the dispersive approach. 
We show that there is no contradiction, our bounds restraining the region for £(w) 
allowed by this latter method. 
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